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Introduction 

Current laser machining using nano-second pulses cannot produce features 
smaller than the spot size due to thermal interactions with the material [1], This paper 
presents arguments and a demonstration that by using femtosecond pulses the spatial 
resolution of the ablation/machining process can be considerably less than the wavelength 
of the laser radiation used to produce it. We have produced laser ablated holes with 
diameters of 10% the laser beam spot size. 



The basic premise for explaining our observations is found in (a) the solid state 
thermal diffusion length for the deposited optical energy, (b) the absorption depth for the 
optical radiation and (c) the time necessary for a vapor to develop at the surface of the 
material. 

For ultrafast pulses in metals the thermal diffusion length, lth=(Dt)^/^ (where D 
is the thermal diffusivity and t the pulse time), is significantly smaller than the absorption 
depth (1/a), where a is the absorption coefficient for the radiation. For our experimental 
conditions with wavelength of 800 nm and 200 fs pulses on gold, the absorption depth is 
275 A with a diffusion length of 50 A. In the case of nanosecond pulses the diffusion 
length, which is on the order of 10 fim, is much longer than the absorption depth, 
resulting in thermal diffusion being the limiting factor in feature size resolution. 
Empirical evidence for the existence of these two regimes is exhibited in Figure 1, Here, 
both experimental [2] and theoretical [3] ablation thresholds are plotted as a function of 
pulse width. An arrow at approximately 7 picoseconds delineates the point at which the 
thermal diffusion length (Uh) is equal to the absorption depth (1/a), It is clear from the 
figure that this is the point at which the ablation threshold transitions from a slowly 
varying or nearly constant value as a function of pulse width to one that is dramatically 
dependent on pulse time. In further suppon of these arguments, it has been demonstrated 
[4] that the electron ihermalization time for laser deposited energy in gold is on the order 
of, or less than, 500 fs and the electron-lattice interaction rime is 1 ps. 

The consequences of this for ultrafast laser pulses is that the energy is contained 
witiiin the beam spot. In fact for energies at or near the threshold for ablation, the spatial 
profile of the laser beam will determine the size and shape of the region being ablated* 

We have performed an experiment to measure the amount of recombination light 
produced as a function of the fluence impinging on a gold film. The technique involved 
is explained in detail in a similar study based on transparent materials [5]. We take the 
intensity of the light to be proportional to the amount of material ablated. In Figure 2, the 
material removed is plotted as a function of fluence. A well defined threshold fluence is 
observed at which material removal is initiated. It is expected that by having only a small 
fi^ction of the gaussian beam where the fluence is greater than the tiireshold, the ablated 
region can be restricted to this small area. In Figure 2, Ra is the radial position on the 
beam where the fluence is at threshold. Ablation, then, occurs only v/ithin a radius Ra- It 
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is evident that by properly choosing the incident fluence, the ablated spot can in principle 
be smaller than the spot size, Rs. Although the data for a 150 fs pulse is shown, this 
threshold behavior is exhibited in a wide range of pulse widths. However, sub spot size 
ablation is not possible in the longer pulse regimes, due to the dominance of thermal 
diffusion. 



Our laser source is a 800 nm TirSapphire oscillator whose pulses are stretched by 
a grating pair, amplified in a regenerative amplifier operating at 1 kHz, and finally 
recompressed by another grating pair. We can obtain pulse widths from 7 ns to 100 fs. 
The beam is focused with a lOx objective, implying a theoretical spot size of 3.0 \im. 
Figure 3 shows an SEM photo-micrograph of ablated holes obtained in a silver film on 
glass, using a pulse width of 200 fs and a pulse energy of 30 nJ (fluence of 0.4 J/ cm^). 
Two holes of diameter approximately 0.3 ^im arc clearly visible. Similar results have 
been obtained in aluminum. 

These results suggest that by, producing a smaller spot size which is a function of 
numerical aperture^and wavelength, even smaller hole s can be machined. Wejiave 
demonstrate the ability to generate the fourth harmonic (200 nm) using a nonlinear 
crystal. Thus by u sing astronger objective lens along with the 200 nm light, holes with 
diameters of 200 angstroms could in pan 

Conclusion 

It has been demonstrated that sub-wavelength holes can be machined into metal 
surfaces using femtosecond laser pulses. The effect is physically understood in terms of 
the thermal diffusion length, over the time period of the pulse deposition, being less than 
the absorption depth of the incident radiation. The interpretation is further based on the 
hole diameter being determined by the lateral gaussian distribution of the pulse in relation 
to the threshold for vaporization and ablation. 
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Pulse width (ps) 
Figure 1 . Ablation threshold vs. pulsewidth in gold. 




Ruence (a.u.) 

Figure 2. Sub spot size machining at ablation threshold in gold. 

Fth is the threshold fluence needed to initiate material removal. 
Rs Is the spot size of the inddent beam and Ra is the radius of the 
ablated hole. 
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Figure 3. SEM view of two 0.3 |im holes machined in silver film. 
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Laser-Induced Breakdown by Impact Ionization! 
in Si02 with Pulse-widths from 7 ns to 150 fs 



D. Du, X. Liu, G. Korn, J. Squier, and G. Mourou 
Center for Ultmfast Optical Science, University of Michigan, Ann Arbor, MI 48109-2099 

Abstract 

Results of laser-induced breakdown experiments in fused silica (Si02) employ- 
ing 150 fs - 7 ns, 780 nm laser pulses are reported. The avalanche ionization 
mechanism is found to dominate over the entire pulse-width range. Fluence 
breakdown threshold does not follow the scaling of Fth ^ when pulses* 
are shorter than 10 ps. The impact ionization coefficient of SiOj is measured 
up to ~ 3 X 10® V/cm, The relative role of photoionization in breakdown for 
ultrashort pulses is discussed. 
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Laser-induced breakdown (LIB) in optically tr«insparent materials has been studied ex- 
tensively since the laser Wcis invented in the 1960's [1]- [4]. LIB mechanisms have been 
investigated for laser pulse-widths down to tens of picoseconds. LIB remains an important 
topic because of its role in a wide range of high-power laser applications, where damage 
to optical components due to LIB often is the ultimate restriction on system performance. 
Increasingly, high laser power is achieved. with moderate amount of energy joules) by re- 
ducing the pulse-width from picoseconds to femtoseconds. Our present work is to investigate 
LIB in this important regime. 

It is a well-established fact that the breakdown threshold depends on the pulse-width 
of the laser pulses. An empirical scaling law of the fluence breakdown threshold for longer 
pulses (Tp > 30 ps) exists [6]: i^th ^/t^. Until recently, the shortest pulses available for LIB 
studies were in the range of tens of picoseconds. However, new techniques in ultrashort laser 
pulse generation and amplification, such as the chirped-pulse amplification (CPA) [7], have 
pushed the available pulse-widths to femtosecond regime. More importantly, CPA readily 
allows one to vary the laser pulse-width continuously — up to hundreds of picoseconds — 
from a single laser system based on this principle. The setup of an LIB experiment using 
such a laser can therefore have a wide range of laser pulse-widths, without changing other 
important parameters, such as the focusing of the pulse onto the s£mfiple. 

We have performed a series of experiments to determine the LIB threshold as a function 
of laser pulse-width between 150 f s - 7 ns, using a CPA laser system. The short-pulse laser 
used in our experiment was a 10-Hz Ti:sapphire oscillator-amplifier system based on the 
CPA technique. The laser pulse was focused by an / = 25 cm lens inside the Si02 sample. 
The Rayleigh length of the focused beam is 2 mm. The focused spot-size was measured 
in-situ by a microscope objective lens. The measured spot-size (FWHM) was 26 /im in a 
gaussian mode. The fused silica samples were made from Corning 7940, with a thickness 
of 0.15 mm. They were optically polished on both sides with a scratch/dig of 20-10. Each 
sample was cleaned by methanol before the experiment. Thin samples were used in order 
to avoid the complications of self-focusing of the laser pulses in the bulk. The Si02 sample 



was mounted on a computer-controlled motorized X-Y translation stage. Each location on 
the sample was illuminated by the laser only once. 

Two diagnostics were used to determine the breakdown threshold Fth- First, the plasma 
emission from the focal region was collected by a lens to a photomultiplier tube with appro- 
priate filters. Second, the change of transmission through the sample was measured with 
an energy meter. Visual inspection was performed to confirm the breakdown at nanosecond 
pulse duration. Fig. 1 shows typical plasma emission and transmitted light signal vs. inci- 
dent laser energy plots, at a laser pulse-width of Tp = 300 fs. It is worth noting that the 
transmission changed slowly at around Fth- This can be explained by the the temporal and 
spatial behavior of the breakdown with ultrashort pulses. Due to the spatial variation of the 
intensity, the breakdown will reach threshold at the center of the^focus, and Becaule of the 
short pulse duration, the generated plasma will stay localized. The decrease in transmitted 
light is due to the reflection, scattering, and absorption by the plasma. By assuming a 
gaussian profile in both time and space for the laser intensity, and further assuming that 
the avalanche takes the entire pulse duration to reach threshold, one can show that the 
transmitted laser energy C/t as a function of the input energy U is given by 

U^^\ (1) 

where k is the linear transmission coefficient. The solid curve in Fig. 1 is plotted using 
Eq. (1), with Uih as a fitting parameter. In contrast, breakdown caused by nanosecond laser 
pulses cuts off the transmitted beam near the peak of the pulses [8], indicating a different 
temporal and spatiaJ behavior. 

Figure 2 shows the fluence breakdown threshold Fth as a function of laser pulse- width. 
From 7 ns to about iO ps, the breakdown threshold follows the scaling law of Fth oc y/r^. 
Breakdown thresholds from earlier work in this pulse width regime are also shown as a 
comparison - it can be seen that our data is consistent with earlier work. However, when 
the pulse-width becomes shorter than a few picoseconds, the threshold starts to increase. 
This increase at shorter pulse-widths may come as a surprise at first, but it still falls into 



the regime of an avalanche dominated process at high field strength, as we will show below. 

The ionization process of a solid illuminated by an intense laser pulse can be described 
by the general equation 

where ne{t) is the free electron (plasma) density, ri{E) is the avalanche coefficient, and E is 
the electric field strength. The second term on the right hand side is the photoionization 
contribution, and the third term is the loss due to electron diffusion, recombination, etc. 
When the pulse-width is in the picosecond regime, the loss of the electron is negligible during 
the duration of the short pulse, 

Photoionization contribution can be estimated by the tunneling rate [9], For short pulses, 
E ~ 10® V/cm, the tunneling rate is estimated to be u; ^ 4 x 10® sec"^, which is small 
compared to that of avalanche, which is derived below. However, photoionization can provide 
the initial electrons needed for the avalanche processes at short pulse-widths. For example, 
our data shows at 1 ps, the rms field threshold is about 5 x 10'' V/cm, The field will reach 
a value of 3.5 x 10^ V/cm (rms) at 0.5 ps before the peak of the pulse, and w 100 sec'^ 
During a ~ 100 fs period the electron density can reach nt[l — exp(-u;A<)] ^ 10^^ 
cm"^, where rit ^ 10^^ is the total initial valence band electron density. 

Neglecting the last two terms in Eq. (2), we then have the case of an electron avalanche 
process, with impact ionization by primary electrons driven by the laser field. The electron 
density is then given by ndt) = no x exp where no is the initial free electron density. 

These initial electrons may be generated through thermal ionization of shallow traps or 
photoionization. When assisted by photoionization at short pulse regime, the breakdown is 
more deterministic, as contrary to that of nanosecond pulses, where the breakdown is more 
statistical [10). According to the commonly accepted condition that breakdown occurs when 
the electron density exceeds nth ^ 10*® cm"^ and an initial density of no ^ 10*^ cm"^, the 
breakdown condition is then given by rjTj, 18 [5], [10]. For our experiment, it is more 
appropriate to use nth 1.6 x 10^^ cm"^, the plasma critical density, hence the threshold 



is reached when riTp 30. There is some arbitrariness in the definition of plasma density 
relating to the breakdown threshold. However, the particular choice of plasma density does 
not change the dependence of threshold as function of pulse duration (the scaling law). 

In our experiment, the applied electric field is on the order of a few tens of MV/cm 
and higher. Under such a high field, the electrons have an average energy of 5 eV, 
and the electron collision time r is less than 0.4 fs for electrons with energy (7 > 5 - 6 
eV [11]. Electrons will make more than one collision during one period of the electric 
oscillation. Hence the electric field is essentially a dc field to those high energy electrons. 
The breakdown field at optical frequencies has been shown to correspond to dc breakdown 
field by the relationship [2] E^{uj) = £??h(l + where lj is the optical frequency 

and r is the collision time. 

In dc breakdown, the ionization rate per unit length, a, is used to describe the avalanche 
process, with tj = a{E)v^ia^ where Vdrift is the drift velocity of electrons. When the elec- 
tric field is as high as a few MV/cm, the drift velocity of free electrons is saturated and 
independent of the laser electric field, Vdrift 2 x 10^ cm/s [12]. 

Thornber [13] has derived an expression for a{E) which is applicable for all electric field 
strengths, which is essential when comparing to our experimental data. The ionization rate 
per unit length of an electron is just cE/Ui times the probability, that the electron has 

an energy > t/,, or a{E) = (eE/Ui)P{E). We denote EkT, Ep, and Ei as threshold fields for 
electrons to overcome the decelerating effects of thermal, phonon, and ionization scattering, 
respectively. When the electric field is negligible, E < EkT, so the distribution is essentially 
thermal, P{E) is simply exp{-Ui/kT). Shockley [14] showed P{E) ~ exp(-const/jE;) for 
EkT < E < Ep. Wolff [15] found P(E) - exp(-const/^2) higher fields {E > Ep). 
Thornber combined the three cases and gave an expression that satisfies both low and high 
field limits: 

This leads to Fth oc JSVp ~ l/xp, i.e., the fluence threshold will increase for ultrashort laser 



pulses when E > y/EpEi is satisfied. 

In Fig. 3, we plot a as a function of the electric field, E. Prom our experimental data, we 
calculate a according to rjTp = 30 and 7? = avdrifi- The solid curve is calculated from Eq. (3), 
using Ei = 30 MV/cm, Ep = 3.2 MV/cm, and EkT = 0.01 MV/cm. These parameters 
are calculated from U = eEl^ where U is the appropriate thermal, phonon, and ionization 
energy, and / is the corresponding energy relaxation length (hr = /p ~ 5 A, the atomic 
spacing, and /»• 30 A, see* [16] and the references there in). It shows the same saturation 
as the experimental data. The dashed line is corrected by a factor of 1.7, which results in an 
excellent fit with the experimental data. This factor of 1.7 is of relatively minor importance, 
as it can be due to a systematic correction, or because breakdown occurred on the surface 
first, which could have a lower threshold. The uncertainty of the saturation value of Vdrift 
also can be a factor. The most important aspect is that the shape (slope) of the curve given 
by Eq. (3) provides excellent agreement with the experimental data. This prompts us to 
conclude that in our experiments, the mechanism of laser-induced breakdown in fused silica, 
using pulses as short as 150 fs and wavelength at 780 nm, is still dominated by the avalanche 
process. 

In summary, we have investigated the pulse- width dependence of the laser-induced optical 
breakdown threshold in fused silica with laser pulse-widths over approximately 5 orders 
of magnitude, from 7 ns to 150 fs. We are able to obtain very good agreement, over the 
entire pulse-width range, on the ionization rate a between experiment and impact ionization 
theory. For picosecond pulses, evidence shows that photoionization is responsible for the 
initial electron generation. The plasma generated during breakdown with ultrashort pulses 
remains localized at the threshold, which could provide precise control of the interaction 
region for material piocessing and medical laser applications when ultrashort pulses are 
used. Our study is also relevant to solid-state microelectronics, as the structures become 
smaller and faster, electron transport properties under high fields and short durations - 
conditions explored in our experiment - will be increasingly important. 
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Figure Captions 



FIG. 1. Scattered plasma emission and transmitted laser pulse as a function of incident laser 
pulse energy. 



FIG. 2i Fluence threshold vs. pulse- width. The .^7]^ scaling holds for pulse- widths down to 10 
ps, shown by the solid line. For longer pulse-widths, previously obtained bulk (squares, Smith [3]) 
and surface (triangles, Stokowski et ai [17]) damage data are also shown. 



FIG. 3. Impact ionization rate per unit distance a determined from experiment and theory. 
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Abstract 

Damage thresholds for corneal tissue are measured as a function of pulse duration, 
revealing strong departure from previous findings. 
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Summary 

We measured damage thresholds for corneal tissue over a range of pulse 
durations (from 150 fs - 7 ns) that could not be examined previously with a single laser 
source, using a Tirsapphire chirped-pulse amplification (CPA) laser [1]. Our results 
depart from previous findings [2.3], which had to use multiple laser sources to investigate 
a similar pulsewidth range. 
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Figure 1. Damage threshold of comea as a function of pulse duration. 
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